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BLOCKAGE CORRECTIONS FOR THREE-DIMENSIONAL-FLOW CLOSED-THROAT WIND 
TUNNELS, WITH CONSIDERATION OF THE EFFECT OF COMPRESSIBILITY 

By John G. Heebiot 


SUMMARY 

Theoretical blockage corrections are presented for a body of 
revolution and for a three-dimensional unswept wing in a cir- 
cular or rectangular wind tunnel. The theory takes account of 
the effects of the wake and of the compressibility of the fluid, 
and is based on the assumption that the dimensions of the model 
are small in comparison with those of the tunnel throat. For- 
mulas are given for correcting a number of the quantities, such 
as dynamic pressure and Mach number, measured in wind- 
tunnel tests. The report presents a summary and unification 
of the existing literature on the subject. 

INTRODUCTION 

When a model is placed in a closed-throat wind tunnel 
there is an effective constriction or blockage of the flow at 
the throat, of the tunnel. The effect of this blockage is to 
increase the velocity of the fluid flowing past the model; 
if the model is not too large relative to the tunnel throat, this 
velocity increment is approximately the same at all points 
of the model so that the model is effectively working in a 
uniform stream of fluid the velocity of which is, however, 
greater than the free-stream velocity observed at some 
distance upstream of the model. It is therefore necessary 
to correct the observed velocity, dynamic pressure, Mach 
number, and other measured quantities for the effect of 
this constriction. This correction is frequently called the 
correction for “solid blockage.” In addition to this solid- 
blockage correction, a correction for “wake blockage” is also 
necessary if the true dynamic pressure and Mach number at 
the model are to be determined. This wake blockage 
arises because the fluid is slowed down in the wake and con- 
sequently must be speeded up outside the wake. A further 
effect of the wake is to produce a pressure gradient which 
must be considered in correcting the drag coefficient. 

Formulas for the solid-blockage correction for a model 
mounted in a two-dimensional-flow wind tunnel are given in 
references 1, 2, 3, and 4. The first-order effects of the com- 
pressibility of the fluid on these corrections are given in 
reference 5, as well as in references 3 and 4. References 
2, 3, and 4 consider also the wake-blockage correction, and 
reference 3 also considers the drag correction due to the 
pressure gradient caused by the wake. The formulas given 
in reference 3 for the solid- and wake-blockage corrections 
will usually be found most convenient whenever the engineer 
is confronted by a practical problem of determining the 
corrections for any configuration met in his experimental 
work. 


The solid-blockage correction for a model mounted in a 
three-dimensional-flow -wind tunnel has been given in a 
number of different forms by different authors. Not only 
do different authors give the correction for the same con- 
figuration in different forms but no one author gives formulas 
which are applicable to both fuselages and wings in tunnels 
of various shapes; for this reason, the engineer confronted 
with a correction problem may have to refer to several reports 
to get the complete solution of his problem. Moreover the 
modifications of the formulas for the first-order effects of 
fluid compressibility are given incorrectly in some cases. 
References 1 and 2 give a formula for the solid-blockage 
correction for a body of revolution in a circular or rectangular 
tunnel for the case of incompressible flow. In references 5 
and 6 the effect of the compressibility of the fluid on this 
correction is discussed, but the result given is incorrect. 
Reference 4 gives a formula for the solid-blockage correction 
for a body of revolution and for a three-dimensional wing in 
a 7- by 10-foot wind tunnel. The modification for com- 
pressibility is correct for the wing but wrong for the body 
of revolution. Reference 7 gives a formula for the solid- 
blockage correction for a body of revolution in a circular 
tunnel, correctly taking account of the effect of the com- 
pressibility of the fluid. References 8, 9, and 10 give a 
formula for the solid-blockage correction for any body in a 
circular wind tunnel together with the appropriate constants 
for a body of revolution and for a rectangular wing having 
various span-to-diameter ratios; the modification of fbis 
formula to take account of compressibility is correctly given. 

It is clear that no one report gives all the necessary 
formulas together with the appropriate constants and com- 
pressibility modifications to enable the engineer to calculate 
the solid-blockage correction for any case with which he 
may be confronted. Moreover, when the results of two or 
more reports overlap, the forms are frequently different so 
that it is not obvious whether the results are in agreement. 
It is the purpose of the present report to summarize and 
extend the results of the previously mentioned reports. 
Formulas are given for the calculation of the solid-blockage 
correction for a body of revolution or a three-dimensional 
unswept wing in a circular or rectangular tunnel. These 
formulas contain tw r o constants, one depending on the shape 
of the body and the other on the shape of the tunnel and 
the ratio of wing span to tunnel breadth. (This ratio may 
be taken to be zero for a body of revolution.) Values of the 
first constant for various bodies of revolution and for a 
number of frequently encountered wing-profile sections are. 
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given. Values of the second constant, for a circular tunnel 
and for rectangular tunnels o_f a number of commonly en- 
countered breadth-to-height ratios are given for various 
wing-span-to-tunnel-breadth ratios. Some of these values 
have been taken from the previously mentioned reports; 
whereas others appear for the first time in the present report. 
The discussion is limited to bodies centrally located in the 
wind tunnel. 

The wake-blockage correction for a model mounted in a 
three-dimensional-flow tunnel is given in references. 4, 8, 9, 
and 10. For the case of incompressible flow the formulas 
are in agreement, but the modification to take account of 
compressibility given in reference. 4 differs from that given 
in references 8, 9, and 10. This matter is discussed in the 
present report. In reference 11 the corrections for the 
pressure gradient due to the wake are given with the correct 
compressibility factors. 

All the final correction formulas together with directions 
for their use are given in the final section entitled “Conclud- 
ing Remarks.” Mathematical symbols are defined as intro- 
duced in the text. For reference, a list of the more important 
symbols and their definitions is given in appendix B. 

SOLID BLOCKAGE IN INCOMPRESSIBLE FLOW 

In studying the flow over a thin airfoil, of small camber at 
a small angle of attack it has been shown that the effects, of 
camber and thickness may be considered independently. In 
treating the problem of wall interference, it is again con- 
venient to consider the thickness and camber effects sep- 
arately. The camber effect, as pointed out in reference 3, 
contributes nothing to the blockage correction. Conse- 
quently it suffices to determine the blockage correction for 
symmetrical bodies at zero .angle of attack. This means 
that for wings it is necessary to consider only the base profile 
of the airfoil, the base profile being defined as the profile the 
airfoil would have if the camber .were removed and the re- 
sulting symmetrical airfoil placed at zero angle of attack; 
bodies of revolution need be considered only at zero angle of 
attack. 

RECTANGULAR TUNNEL 

Three-dimensional wing. — The blockage correction for a 
three-dimensional wing in a rectangular tunnel is considered 
in reference 4. Numerical values are given only for a wing 
of 6-foot span in a 7- by 10-foot wind tunnel. The method 
may, however, be applied to any rectangular tunnel and any 
span-to-breadth ratio. 

Consider a rectangular tunnel of height H and breadth R. 
Suppose that the wing span 2s is in the R-direction so that 
2 s/B is the span-to-breadth ratio. As in reference 4, let 
the wing be represented by a series of finite lines of sources 
and sinks. (See fig. 1.) The strengths of these sources 
and sinks are assumed to depend on tire airfoil profile and 
the determination of these strengths, which is a two- 
dimensional problem, is explained presently. If each line 
retains the same strength from one end to the other, the 
wing section cannot be exactly constant. It will thin 
down at the extreme tip and the plan form will not be 
exactly rectangular, but neither of these features is such as 
to detract from its usefulness for the present purpose, which 



Fiogke 1.— Image system for three-dimensional wing In rectangular tunnel. 


is to represent an actual wing sufficiently well to enable a 
calculation to be made of the velocity along the tunnel 
axis “induced” by images of the wing. It is this velocity 
induced by the images which represents the effect of the 
tunnel walls on the velocity at. the model, and which is the 
solid-blockage correction. 

Consider the image line source CD of strength Q per unit 
length. The velocity potential of a three-dimensional 
source oi strength Q5y (volume per unit time) is —Q5y[4.wr. 
It follows that the component velocity along the tunnel 
center line induced at A by the source element QSy is 


5 u 



ty 


where g and r are defined in figure 1. Putting r*= 
(mB—y) t +n 2 H 2 +g t and integrating from —s to 1 .gives 

a u „ ■ Qg- f mB+ * 

d 4 -k{ti 2 IP + <7 S ) [_ yh 2 // 2 + g % -1- (mB s) ! 

1 (d 


for the single line source. CP. Now CD is one image of one 
finite line source used to represent the wing in the tunnel. 
In order to find the velocity induced by all the images of this 
particular line source, it is necessary to add the results 
obtained from equation (1) by giving m and n all positive 
and negative integral values except m=n= 0. Thus for a 
single line source. 

v. —Q& T'' J- P mB-^s 

rctP+g* 1 _ JrfTP+g'+imB+sy 


mB—s 

+ g* + (mB—sy 


i 


( 2 ) 
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where the prime denotes that the term m=n=0 is to be 
omitted from the summation. Equation (2) may be re- 
written 


where 


UA = 2^P-r 



l 9 $ R \ '•v-V/ 1 

a \H’ B’ H)~2s ^ n 2 +(g[H) 2 

m-rsIB 




i +(g/Hf+(m+sfBf(BjH) 
m—s/B 


^+(glH) t +(m -simBm 


2 


(3) 


where 

f maximum thickness of airfoil 
C cross sectional area of tunnel 

If T denotes the volume of the wingj then V=2sctxi where 
K t < 1 depends on the shape of the base profile. Equation (7) 
may be rewritten 


or 


where 


A JT_K 1 tV_ 

' JJf ~ <03/2 


AJJ' K 2 r2sct 
~U r ~ C Vi 


( 9 ) 



( 10 ) 


It is convenient to define 


T 


( 


g_ ± B\_ 1 ( B V /2 ( g $_ B\ 
H! B’ h )~ 2 \tH) ff \H’ B' HJ 


so that 


,, Vss 
A- f BH ?' 2 


1 B\ 
\H’ B’ HJ 


(4) 


Now a and r depend only very slightly on g[H and so it 
usually suffices to take g/H= 0 in the evaluation of these 
quantities. Any wing profile can be represented by a suitable 
distribution of sources and sinks along the chord. It follows 
that the total induced velocity due to the wing images is 
obtained by summing over this distribution and is approxi- 
mately 


\ Tjf s ~Qg 

lL ~{BH) Vi 



( 5 ) 


It may be noted here that setting g[H= 0 in the evaluation 
of u and r is equivalent to representing the wing by a line 
doublet of strength ~ZQg (analogous to references 1 and 2) 
instead of by a distribution of sources and sinks. 

The quantity 2<2p of equation (5) can be determined 
approximately from the wing profile and this determination 
is a two-dimensional problem. From reference 2, but with 
the notation of reference 3, there is obtained 


3Gf-fAe*Z7' (6) 

where 

c airfoil chord 

TJ' apparent free-stream velocity at airfoil as determined 
from measurements taken at a point far ahead of 
model 

A a factor dependent on shape of base profile 


Substitution of equation (6) into equation (5) yields 

AiL r/ 2sc 2 w 3/2 , /„ s 
U' ~(BH) m 16 V’ B’H/ 

_2sctr 3li A / s B\ 

~Cv* 16 t/c T V U ' B’ HJ 


(7) 


A1 

1 16 t/c k, 


A 

1 16 t/c 


(ID 

( 12 ) 


It is clear that r depends only on the timuel shape and the 
wing-span-to-tunnel-breadth ratio; whereas K lt K* depend 
only on the shape of the base profile. 

The factor A can be determined for any base profile from 
the relation (references 2 and 3) 

A== vl 1 V i ^ Vl +(.dy t /dx)*d (f) (13) 

where 

y, ordinate of base profile at chord wise station x 
dyjdx slope of surface of base profile at x 
P base-profile pressure coefficient at z in an incompres- 
sible flow 

Values of A for a number of base profiles are given in refer- 
ence 3. Thus the value of K 2 can be calculated from equa- 
tion (12). The value of is immediately found from the 
area of the base profile which may be calculated, for example, 
by a numerical integration from the ordinates of the base 
profile. As soon as K t is known, Ki can then be calculated 
from equation (11). The evaluation of r requires the sum- 
mation of the infinite series in equation (3). The summation 
of this series is explained in Appendix A. 

Values of r for rectangular tunnels of various breadth-to- 
height ratios and for various wing-span-to-tunnel-breadth 
ratios are given in table I and figure 2. Values of K t and K 3 
for various base profiles are given in tables II and III. The 
choice between the two formulas (S) and (9) is entirely a 
matter of convenience and should be decided in the light of 
the available data. 

Body of revolution. — The blockage correction for a body 
of revolution in a rectangular t unn el is considered in refer- 
ences 1, 2, and 4. In reference 4 numerical values are given 
for some average streamline body of revolution in a 7- by 
10-foot wind tunnel; whereas in references 1 and 2 numerical 
values are given for prolate spheroids and Rankine Ovoids 
in square and duplex (breadth equal to twice the height) 
wind tunnels. Either the method of reference 4 in which 
the body of revolution is represented by a suitable distribu- 
tion of sources and sinks along its chord or the method of 
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TABLE I— VALUES OE r FOR VARIOUS TUNNEL SHAPES 
AND CONFIGURATIONS 


Tunnel shape ^ N \ 

Body of 
revolu- 
tion 

3-dimensional wing with spa ri-to -breadth ratio 

f- 


I- 0 - 50 ' 


§=ux> 

Circular 

Square 

/B/ 17=10/7 
BIH= 7/4 
BIH- 2 

Rectangular. ..j B/H-712 
)BI 11=2/7 

73/11=1/2 

73/17=6/10 

lB/U=7/10 

0.797 

.812 

.863 

.946 

1.028 

1.729 

1.729 

1.028 

.922 

.863 

0.797 
.812 - 
.863 
.946 
1.028 
1.729 
1.729 
1.028 
'..922 
.863 

0.812 

.818 

,864 

.941 

1.017 

1.630 

1.783 

“’932 

0.828 

.836 

.866 

.930 

.990 

1.436 

1.896 

.977 

0.859 

.874 

.884 

.923 

.967 

1.204 

2.196 

1.063 

0.951 
.916 
.937 
.9G2 
1.160 
2.665 
1.434 
1.230 
1. 110 


TABLE II.— VALUES OF Ki FOR VARIOUS BASE PROFILES 


\ 

fo.xx\ 

Ellipse 

Joukow- 

ski 

section 

Conven- 

tional 

NACA 

sections 

OOXX 

NACA low-drag sections 

I6-0XX 

64-OXX 

6HXX 

66-0XX 

0.06 

0,938 


0.941 

0.938 

0.962 

0.965 

0.987 

.09 

.965 

6.991 

.972 

.962 

.961 

.967 

.984 

.12 

.993 

1.016 

1.005 

.997 

.987 

.989 

1.006 

.15 

1.019 

1.045 

1,035 

1.028 

1.019 

1.020 

1.032 

.18 

1.046 

1.068 

1.063 

1. 062 

L 047 

1.051 

1.057 

.21 

1.072 

1.083 

1.090 

L 090 

1.073 

1.079 

1.079 

.25 

1. 108 


1.128 

L 128 




.30 

1. 152 







.35 

1. 196 






» 

.50 

1.329 


. 

" " 




1.00 

1.772 





, 











TABLE III— VALUES OFF, FOR VARIOUS BASE PROFILES 


-=o.xx\ 

« \ 

Ran- 

kine 

Oval 

Ellipse 

Joukow- 
skl sec- 
tion 

Conven- 

tional 

NACA 

sections 

OOXX 

NACA low-drag sections 

■ 

16-OXX 

04-0XX 

65-0XX 

66-OXX 

0.06 


0.737 


0.644 

0.690 

0.615 

0.632 

0.679 

.09 

0.913 

.758 

0.59ft 

.665 

.708 

.611 

.630 

.673 

.12 

.928 

.780 

.615 

.687 

.734 

.624 

.641 

.684 

.15 

.935 

.800 

.633 

.708 

.756 

.640 

.657 

.698 

.18 

.953 

.822 

.652 

.727 

.781 

.654 

.673 

.712 

.21 

.961 

.842 

.670 

.746 

.802 

.666 

.686 

.723 

.25 

.979 

.870 

.682 

.771 

.830 

.678 

.697 

.732 

.30 

1.002 

.905 

.726 





p 


.36 

1.043 

.940 

.763 






.50 

1.176 

1.044 

.876 






1.00 

1.392 

1.392 











- 





references 1 and 2 in which the body is represented by a 
doublet of suitable strength at its center may be used to 
obtain results for any streamline body of revolution in any 
rectangular tunnel. Since the method of reference 4 was 
used for the three-dimensional wing, it is instructive to use 
the doublet method of references 1 and 2 for the body of 
revolution case, although both methods give the same results. 

Consider a body of revolution of maximum thickness t and 
length c centrally located in. a rectangular tunnel of breadth 
B and height II. As in reference 2 the body may be repre- 
sented by a doublet of strength n given by the equation 

\W (14) 
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(a) Variation of t with ratio of wing span to tunnel breadth, for various tunnel sham-s. 


Figure 2.— Values of r for various tunnel shapes and configurations. 



Fioure 2.— Concluded. 


the summation being taken over all positive and negative 
integral values of m and n except m=n— 0. If gfll—O 
equation (3) may be rewritten 

s B 1 r(m+s/B)F mn —(m—8jH)E mn ~] 

'V’B’HrTs^ izzp- J 

-sr-v/ ‘Inn 

2 - J E mn F mn [(m +s/B)F mn +(m - s/B) E m ,\ 

( 10 ) 

where the quantities E mn and F mn , which are introduced for 
convenience, are defined by the equations 

F mn = Jn>+(m- S /B)KB/Hy 

It follows that 


where X is a constant depending only on the shape and fine- 
ness ratio of the body. The velocity induced at the model 
by the tunnel walls is the same as that induced by a doubly 
infinite array of images of the doublet and is given by 


<7 


(o,o,i)- S ' 


1 


( 17 ) 


Substitution of equation (17) into equation (15) yields 



1 

(n 2 H 2 +m 2 B ^ 2 


AiZ7'= 


4 xfP 



(15) 
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If jit is replaced by its value from equation (14) and equation 
(4) is used, there is obtained 


A JJ' *«* 

U' 8 (BH) Slt T 



B\ ct 2 r 3 ' 2 \t 
’ H)~C W 8 c T 


(Oi 0) |) (18) 


If T’ denotes the volume of the body of revolution, then 
I’= K 2 ct 2 where k 2 < v/4 depends on the shape of the meridian 
section of the body. (For a right circular cylinder whose 
meridian section is clearly a rectangle Ki=7r/4.) Equation 
(18) may be rewritten 



AJJ' K 3 tV 
U' ~ C 3 '* 

(19) 

or 

where 

A JJ' K 4 rcf- 

JJt £*/2 

(20) 


t=t( o, o, J) 

(21) 


sr 3/2 Xf 1 
8 c Ki 

(22) 


„ , r 3/2 Xt 

Ki 8 c 

(23) 


TABLE IV —VALUES OF K t FOR VARIOUS BODIES OF 
REVOLUTION 


\ 

t \. 
C 

Prolate 

spheroid 

Rankine 

Ovoid 

NACA 

111 

NACA 

133 

Fuhr- 

TORxm 

source- 

sink 

body 

Sym- 

metric 

Fuhr- 

mitnn 

source- 

sink 

body 


0.08 


0.913 

0-902 





.10 


.923 

.909 



_ 


.12 


.932 

.916 


HI 



.14 

.917 

.941 

.924 





.16 

.924 

. 949 

.931 




■■ 

.18 

.931 

.957 

.939 


0.933 

0.032 

0.925 

.20 

.938 

.964 

.945 




HI' 

.24 

.954 

.979 

.964 





.30 

ill . il 

1.000 



0. 983 




.40 

1.025 








.50 

1.072 






HI 

1.00 

1.329 


— 

— 
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TABLE V— VALUES OF K t FOR VARIOUS BODIES OF 
REVOLUTION 


£ \ 
c \ 

Prolate 

spheroid 

Rankine 

Ovoid 

NACA 

111 

NACA 

133 

Fuhr- 

maim 

source- 

sink 

body 

Sym- 

metric 

Puhr- 

mann 

source- 

sink 

body 

as, 

-i,o) 

body of 
refer- 
ence 14 

0.0S 

0.472 

0.670 

0. 405 





.10 

.474 

.668 

.409 




mimm 

.12 

.477 

.665 

.413 




&b|^h| 

.14 

.480 

.662 

.417 




hshi 

.16 

.484 

.659 

.422 




mmm 

.18 

.438 

.655 

.428 


0.417 

0.414 

0.436 | 

.20 

.491 

.650 

.433 


.24 


.640 

.445 



9HBHH 

■|HH 

.30 

.514 

.624 

0.462 



S535p5JB 

.40 

.537 

.607 





.50 

.562 

.609 





BHHjB 

1.00 

.696 

.696 












■1 


It should be noted that r for the case of the body of revo- 
lution is the same as for the limiting case of a wing when the 
span approaches zero. As pointed out in reference 2, X may 
be calculated for any body of revolution whose pressure dis- 
tribution is known. The necessary formula is 

x=4 (f)T 0)’-^ (f) VV 


where 

y radius of body at ckordwise station x 
dy[dx slope of meridian section at x 
P pressure coefficient at x in incompressible flow 
Values of X for prolate spheroids and Rankine Ovoids are 
given in references 1 and 2. As soon as X is known for a 
body, Ki can be calculated at once from equation (23). The 
value of k 2 may be found from the volume of the body of 
revolution which may be calculated for example by a numer- 



NACA Hi; tic = 0.20 



Prolate spheroid; t/c - 0.20 


Rankine ovoid; t/c • CXSO 




NACA til; t/c -0.18 



Fuhrmann source-sink body; t/c =0.10 



Symmetric Fuhrmann source-sink body ; t/c - 0.1 8 



(3,6,-/,OJ Body of reference 14% -t/c =G.1S 
Figure 3.— Sample bodies of revolution. 
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ical integration, and K s can be. calculated from equation (22) 
as soon as k 2 is known. 

Values of r for rectangular tunnels of various breadth-to- 
lieiglit ratios are given in table I and figure 2. Values of 
Kz and for various bodies of revolution are given in 
tables IV and V. Some of these bodies are drawn in figure 
3. Unfortunately these tables are rather incomplete and 
moreover fuselage shapes used in practice are extreme^ 
varied. However, in table IV it is observed that the values 
of Kz do not depend very strongly on the shape of the body, 
although they do depend on the thickness ratio. For this 
reason it appears that for most fuselages it will be sufficiently 
accurate to use the values of Kz given for the NACA 111 
bodies. As the values of K* giVen in table V are more de- 
pendent on the body shape, it is recommended that equation 
(19) and table IV be used in preference to equation (20) and 
table V whenever possible. ... 

CIRCULAR TUNNEL 

Body of revolution. — It is convenient to consider the case 
of a body of revolution before considering the case of the 
three-dimensional whig because more attention has been 
given to the former by other authors. _ (See references 1, 2, 
7, 8, 9, and 10.) It will now be shown that the blockage 
correction is again given by equations (19) and (20) where 
r has a value appropriate to a circular tunnel. Since K» and 
Ki depend on the model and not on the tunnel, they are 
still given by equations (22) and (23). 

The most convenient starting point is the formula of 
references 1 and 2, namely, 

< 25 ) 

where S m is the maximum cross-sectional area of the model, 
It is only necessary to note that 


where 

Xv factor, depending on model shape 
tv factor depending on tunnel shape 
V volume of model 
B diameter of wind tunnel 

It is of interest to show that equations (19) and. (20) can be 
deduced also from equation (27) showing the latter to be 
equivalent to equation (25). From reference 8, there is 
obtained 



4 

Tir=— r 
■K 


Since also C=rB 2 / 4 equation (27) yields at once 

A zU' 4 \ T7 /V /S 1\ Xfx^cf 2 

W V4 1 A// V8 C z > i J r c 8 C 3 '* 


This is the same as equation (26) which yields equations 
(19) and (20) directly. 

Three-dimensional wing. — The blockage correction for a 
three-dimensional wing in a circular tunnel is given in refer- 
ences 8, 9, and 10, the formula being of the same form as for 
a body of revolution in a circular tunnel, namely, equation 
(27) where r v and Xv have values appropriate to the three- 
dimensional wing. From reference 8 (using the notation of 
reference 3 instead of reference 2) there is obtained 


Xv 


TT 2 sc 2 

8 r 


A 


(28) 


4 




T 


(29) 


Since also C=Bir ! /4 equation (27) yields at once 


S m =T?/ 4 
V=K t ct 2 


It follows at once that . 

A iU' , f 3 Xt 77 3 ' 2 d 2 

IT ~~ tX 8 c 8 W 2 


(26) 


Substitution from equations (22) and (23) reduces this to 
equations (19) and (20), respectively. It should be noted 
that r of references 1 and 2 is identical with r of the present 
report. 

The value of r is given in table I and figure 2. Values of 
Kz and K t are given in tables IV and V. Again equation 
(19) is preferable to equation (20). 

The result of reference 7 fails to take into consideration, 
the body shape and so it is not very useful except for less 
exact calculations,. The result of references 8, 9, and 10 is 
presented in a different form, namely, 


Aii7' Z 

jjr — AvTv 


(27) 


A 1 I^_/t2££ 2 \/4 1 V 2 td A 

U’ \8 V )\* T ) \ 8 C ii2 )~C in t/c 16 r 

This is the same as equation (7) which leads directly to 
equations (8) and (9) with the same definitions of Kz ami 
Kz given in equations (11) and (12); r is, however, given by 
equation (29) where t v is obtained from references 8, 9, or 10. 

Thus equations (8) and (9) may also be used for circular 
tunnels provided only that the appropriate values of r are 
used. Values of r are given in table I and figure 2. Values 
of Kz and K-i are given in tables II and III. 

SOLID BLOCKAGE IN COMPRESSIBLE FLOW 

In the preceding section the solid-blockage corrections 
have been determined under the assumption that the fluid 
is incompressible. It is now necessary to determine the 
modifications required in these formulas to lake account of 
the effects of the compressibility of the fluid. The methods 
of references 12 and 13 are very convenient for this purpose. 
As the required modifications are given iucorrectly in refer- 
ences 5 and 6, partially incorrectly in reference 4, and 
correctly in references 7, 8, 9, and 10, it appears worth while 
to give some discussion of the matter. 
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For the purpose of deducing the properties of a compressi- 
ble flow from those of a corresponding incompressible flow 
the so-called “Extension of the Prandtl Rule,” which was 
first given in reference 12 and repeated as Method IV in 
reference 13, is probably of most general application; but 
the other methods of reference 13 are sometimes more con- 
venient for certain problems. The Extension of the Prandtl 
Rule may be expressed in the following manner: 

The streamline pattern of a compressible flow to be cal- 
culated can be compared with the streamline pattern of an 
incompressible flow which results from the contraction of the 
y and z axes including the profile contour by the factor 
-\j 1 — d/ 2 (d/=free-stream Mach number) (x axis in the 
direction of the free stream). In the compressible flow the 
pressure coefficient as well as the increase in the longitudinal 
velocity are greater in the ratio 1/(1— JP) and the stream- 
line slopes greater in the ratio 1/^1— iP than those at the 
corresponding points of the equivalent incompressible flow. 

Since formulas (8) and (19) for the three-dimensional wing 
and the body of revolution have the same form and also 
apply to both rectangular and circular tunnels, it suffices to 
determine the modification due to compressibility for them. 
Let the subscript c refer to compressible flow and the sub- 
script i to the corresponding incompressible flow. If Y c is 
the volume of the model in the compressible flow, then V<, 
the volume of the model in the corresponding incompressible 
flow, is given by 

V t =[l-Qir]Vc 

where M' is the apparent free-stream Mach number at the 
model as determined from measurements taken at a point 
far ahead of the model. Also if C e is the cross-sectional 
area of the tunnel in the compressible flow, then Ci the 
cross-sectional area of the tunnel in the incompressible flow 
is given by 

c t =[i-{M r nc c 

t is unaffected by the transformation and the effect on Ki 
and K 3 is sufficiently small that it may be neglected. From 
the Extension of the Prandtl Rule it follows that 


/A t C7'\ _ 

- 1 , 


V u r K 

1 — LI/') 2 

K U' A 


1 K,t\\ 1 KjtV c 

1— LI/') 2 C?'* [1-Q/') 2 ] 3 ' 2 C c vl 

where j denotes the numbers 1 or 3. Since equations (8) 
and (9) are equivalent as are also equations (19) and (20), 
it follows that in all cases it is only necessary to multiply 
the blockage corrections given by these formulas by 
[1 — (M') 2 ] _3/2 in order to take account of the compressibility 
of the fluid. 

As a check it is useful to determine the compressibility 
modification for the case of a body of revolution in a circular 
tunnel by Method II of reference 13, since the derivation 
is so simple by this method. Both the body shape and 
the longitudinal velocities are the same in the corresponding 
compressible and incompressible flows; only the tunnel 

dimensions are altered by the factor -\/l — (.!/') 2 so that 
(?i=[l— (M')*](7 e . There is obtained 


/A,Z7'\ K % rV, 1 K z tV c 

\ U' ) e \ U' Jt C 3 ' 2 [l-p/OT' 1 C: 3 ' 2 

as before. 

WAKE BLOCKAGE IN COMPRESSIBLE FLOW 

The blockage due to the wake of a model in a two-dimen- 
sional-flow tunnel is discussed in some detail in reference 3. 
Much of this discussion is applicable without change to the 
case of a three-dimensional- flow tunnel. The fundamental 
idea of replacing the model and its wake by a source (in tills 
case a three-dimensional source) of suitable strength located 
at the position of the model can be used again and in fact the 
determination of the source strength Q can be carried out in 
exactly the same manner. The result is identical with that 
of reference 3, namely, 

q= p' u 'CYS [1 + (y _ 1)( u , )2] (3 0) 


where 

Q mass flow of source rather than volume flow as used 
previously 

p' mass density of fluid at point far upstream 
CV uncorrected drag coefficient referred to apparent dy- 
namic pressure q' 

S area on which drag coefficient is based 
7 ratio of specific heat of gas at constant pressure to 


specific heat at constant volume 



In equation (30) powers of M' higher than (J/') 2 have been 
neglected. 

Consider now a rectagular tunnel of height H and breadth 
B. The tunnel walls are replaced by a doubly infinite array 
of sources of strength Q at distances mB to the side and nH 
above and below the position of the model. 

By making use of Method II of reference 13, it is readily 
shown that a three-dimensional source of strength Q (mass 
per unit time) in a uniform flow of compressible fluid will 
induce at the point the coordinates of which are x, y, z rela- 
tive to the source a streamwise velocity. 


Aw 


Q x 

4^ [z’-Ki-AW+x 2 )] 3 ' 2 


where the uniform flow is in the x direction and p and M are 
the density and Mach number, respectively, of the undis- 
turbed stream. 

It follows that the streamwise velocity A 2 I7' induced at a 
point on the center line of the tunnel by the entire system 
of images is 

AiU ' = I^y ^ {x 2 + [1 - (A/') 2 ] (m 3 £ 2 + n 2 # 2 )} 3 ' 2 


where p' and M' are the density and Mach number of the 
undisturbed flow in the tunnel and the summation is taken 
over all positive and negative integral values of m and n 
except 7i=0. The velocity induced by the image 
sources at an infinite distance upstream is 

(A,B T 0-«==Jim {x 2 +[l-(M') 2 ](77i 2 R 2 +7i 2 f? 2 )} s ^ 
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But- conditions far upstream must remain unchanged and to 
achieve this it is necessary to counterbalance this velocity 
by superposition of a uniform flow of. equal magnitude but 
opposite sign. The addition, of this flow at all points in the 
field will result in a speeding up of the general flow at the 
position of the airfoil by the amount 

x 

7 S' {x 2 +[1 _ } 3/2 (3 1} 

The summation of this series can be obtained by the following 
artifice. Substitution of . equation (30) into equation (31) 
and setting M'= 0 gives 

r< / orp r . .... 

(& s U') t =lim ° 8v jy {a .* +ni 2£2 +7l 2 H * )3 /2 (32) 


WAKE PRESSURE GRADIENT IN COMPRESSIBLE FLOW 

The effect of the pressure gradient caused by the wake is 
discussed in reference 3 for compressible flow in a two- 
dimensional wind. tunnel and in reference 1 1 for compressible 
flow in a circular wind tunnel. The method of reference 1 1 
is equally applicable to the. case of a rectangular wind tunnel 
if the appropriate value of r (table I) is used. 

The longitudinal velocity increment due to the effect of 
the tunnel boundaries on the source used to simulate the 
wake has an approximately linear gradient in the stream 
direction at the model location. This linear gradient in the 
velocity is equivalent to a linear gradient in the pressure ns 
is easily seen from the approximate relation. 


But according to references 4, 8, 9, and 10 

I A Z U'\ 1 C D 'S 
V. U' )<~4 BH 


(33) 


In reference 11 it is shown that the gradient in A 5 f r ' for a 
source in a wind tunnel is identically equal to the value of 
A t U' for a doublet hi a wind tunnel. In the notation of 
the present report there is obtained 


Comparison of equations (32) and (33) shows that 


lira S' 


X 2 T 

(z 2 + mW+nWy^BH 


(34) 



For compressible flow this becomes 


If in equation (34) B and II are replaced by Vl— (M') 2 B 
and ^\~(M'fII, respectively, tlierelS obtained 

.. x 2ir 

{ x 2 + [ 1 - {My} } 3/2 [l-(M'¥]BH 

(35) 

Substitution of equations (35) and (30) into equation (31) 
yields 


A Z U' 1-Hr-lKM') 2 C D 'S 1+0.4 (iWyC D 'S 
U' ~ l -{M'f 4BII~ l— (My 4:C { ) 

on setting y=1.4. The preceding discussion is for a rectan- 
gular tunnel, but in references 8, 9, and 10 the same formula 
is given for a tunnel of any shape for the case M'= 0. Con- 
sequently, equation (36) may be taken to hold generally for 
a tunnel of any shape. 

In reference 4 the wake-blockage correction is given cor- 
rectly for incompressible flow but the modification to take 
account of compressibility is given incorrectly, as it is deter- 
mined from reference 5 which, as pointed out in reference 13, 
is incorrect. The effect of compi'essibility on the wake- 
blockage correction is determined in references 8, 9, and 10 
by means of the correct form of .the Extension of the Prandtl 
Rule but there is some question whether this rule is applica- 
ble to this problem. It appears that in using this method 
to determine the effect of compressibility on the blockage 
correction, the effect of compressibility on the source strength 
Q given by equation (30) is overlooked. This explains the 
discrepancy between equation (36) and the formula of refer- 
ences 8, 9, and 10. 


d£ 

dx_ 


i+(7-i my 
2 [1-(M') 2 ] S ' 2 



CVS 

QUt 


1 + 0.4(M')* R Co'S 

“ 2 [1— (M') 2 l J/! V 4 7 W 


The increase in the drag resulting from this pressure 
gradient is equal to the product of the pressure gradient 
by the sum of the actual model volume and the virtual 
volume (reference 2). It should be noted that t he con stants 
Ki and K z of tables II and IV are equal to V(jr/4) limes 
the ratio of the sum of the actual model volume anti the 
virtual volume to the actual model volume, these quantities 
being calculated under the assumption of incompressible 
flow. Tims the increase in drag coefficient caused by the 
pressure gradient is given by 


A On' 


for the wing, and 


AC\' = 


1 + 0.4(M') 2 KxrVmCo' 
[1 — (A/') 2 ] 3/ * W 1 


i+0A(My K t T\\n 0 ' 
[1 -(MW C % > % 


(37) 

(38) 


for the body of revolution. 

It is pointed out in reference 11 that the virtual volume is 
altered by the compressibility of the fluid. Thus equations 
(37). and (38) can be slightly improved if K\ and K z appearing 
therein are corrected to take account of this effect. The 
necessary modification is made by replacing K\ and K z in 
these equations by Iu„ and K Zv , respectively, where 

j ^ = Vf [} +mn (S Ki ~ O]—. 3 0 °> 
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Fifli'R* Linear compressibility correction factor for virtual volume. From reference 11. 

where h (M r ) is given by figure 4 which is reproduced from 
reference 11. It should be noted that the improvement in 
modifying K and K 3 for compressibility will be small es- 
pecially in the case of K% for the body of revolution. It is 
important only for quite high Mach numbers. This addi- 
tional refinement was not made in reference 3. 

It appears that a similar compressibility modification of 
Ki and K 3 in the formulas for the solid blockage should be 


mgde, but the exact modification required is not known. 
However, it is believed to be small. 

CORRECTION OF MEASURED QUANTITIES 

The true velocity U at a model consisting of a body of 
revolution and wing can be obtained from the apparent 
velocity V' by applying the solid-blockage and wake- 
blockage corrections. The true velocity may be written in 
the form 

U=U'(1+K) (40) 

where 

K=K a +K>+K ut (41) 


In equation (3S) K a is the solid-blockage correction due to 
the wing and is given by 


K u 


1 

[1 -(MW 3 


kvr, 

C vi 


(42) 


or 


K a = 


K 2 t2scJw 

[i-ao 2 ] 372 c z!i 


(43) 


Kt is the solid-blockage correction due to 
revolution, being given by 

Tr 1 K 3 t V b 

[i ( l /') 2 ] 3 ' 2 Q V1 


the body of 


(44) 


or 


Kt 


1 K rctft 2 

[l-(!f') 2 ] 3/i C 311 


( 45 ) 


K at is the wake-blockage correction and is given by 


r, 1+0 A(M'?C d 'S 

1 -(!/')* 4 C 


(46) 


It is evident that a correction to the apparent velocity in 
a compressible Sow implies corrections also to the apparent 
density, dynamic pressure, Reynolds number, and Mach 


TABLE VI.— COMPRESSIBILITY FACTORS FOR CORRECTION EQUATIONS 
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number. These corrections are readily obtained on the 
basis of the usual assumption that the flow is adiabatic... It 
is assumed that the correction terms are small compared 
with unity, so that squares and products of these terms may 
be neglected. The analysis follows the lines of reference 3, 
and it is not necessary to repeat the details here. The fol- 
lowing equations are obtained: 

p=p'[l-(I/') 2 Aj (47) 


2 = 2 '{l+[2-(M0 2 ]tf} . (48) 

— O.TCM') 2 ]#} (49) 


M=M'{ 1 + [1+0.2 {M'f]K} (50) 

The drag coefficient must be corrected for the. effect of. the 
pressure gradient due to the wake as well as to refer it to the 
correct dynamic pressure. There is thus obtained 

C»=C D '{ l-[2-(My]K} -ACV-A<V (51) 

where A C D J and A C Dt ' are given by equations (37) and (38).. 

Numerical values of the functions of M' which appear in 
these equations are given in table VI. 

CONCLUDING REMARKS 


Data obtained from tests of three-dimensional models, 
which are small relative to the wind-tunnel dimensions, can 
be corrected for solid and wake blockage and for the pressure 
gradient due to the wake by means of the following equations : 


U=U'(l+K) (40) 

2 = 2'{1+L2 -(MW} (48) 

A=A'{1 + [1-0.7(M') 2 ]A} (49) 

M=M'{l + [l + 0.2(M')*\K} (50) 


C D = Co' { 1 - [2 - (MO 2 }K} -AC d J-AC D i ; (5 1) 

In the preceding equations K is obtained from the following: 


where 


or 


K=K a +K,>+K ut 


. 1 K\t V w 

m ~ [1-(M') 2 ] 8/2 C 3/2 

- 1 7^2 T 2 SCmt y, 

'“ = [l_(M , ) 2 l a/2 


(41) 

(42) 

(43) 


Jr 1 AV \ \ 

[T (J/')2]3/2 QZ/i 

(44) 

r, ’ 1 KiTCbf* 

b ~ [1 - (J/'j+A 

(45) 

„ 1 + 0.4(A/T Co'S 

^ vk ~ 1 -(My 4 c 

(4G) 

a /■> , l + OitM'f/frr/.' 
[1 — (AT) 2 ] 8 '* C 3/1 

(37) 

a /— ( , l + 0.4(Af')*A S TF>Cy 

[l— (M') 2 ] 3/2 c m 

(38) 


In these equations r is a factor which for a body of revolu- 
tion depends only on the shape of the tunnel; whereas for a 
three-dimensional wing r depends on the ratio of the wing 
span to the tunnel breadth as well as on the tunnel shape. 
(The wing span is in. the direction of the tunnel breadth.) 
Values of r are given in table I and figure 2. The constants 
Ki and K 2 for the three-dimensional wing depend only on 
the wing base profile shape and can be calculated by means 
of equations (11), (12), and (13). Values of Ki ami Kt for 
a number of wing profiles are given in tables II and III. 
The constants K s and Jv 4 for the body of revolution depend 
only on the shape.of the body and can be calculated by means 
of equations (22), (23), and (24). Values of these constants 
for a number of body shapes are given in tables IV and V; 
some of these shapes are drawn in figure 3. Since these 
tables are incomplete and fuselage forms are not standard- 
ized as are wing sections, it is recommended that the values 
of Kz given in table IV for the NACA 111 series of shapes 
be used for any fuselage shape which does not differ loo 
greatly from an NACA 111 shape. This implies that equa- 
tion (19) and table IV should be used in preference to equa- 
tion (20) and table V whenever possible. Numerical values 
of the functions of M' which appear in the correction equa- 
tions are given in table VI. 

The constants Ki and K 3 appearing in equations (37) and 
(38) may be modified for compressibility by means of figure 
4. The modified values of Ki and K s are to be used in 
equations (37) and (38) only and not in equations (42) and 
(44). --- 
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APPENDIX A 

SUMMATION OF THE INFINITE SERIES FOR <7f0. a/B, B/H } 


From equations (3) and (16) it is seen that <r(0, s[B, BfH) 
may be written in the alternative forms. 


<r 



1 / m+s/B 
n 2 \ E mn 


m—s/B\ 

Epin / 


(Al) 


or 


where the summation is taken for all positive and negative 
integral values of rn and n except m—n = 0 and the quantities 
E mn and F mn , which are introduced for convenience, are 
defined by the equations 

4»= YV+Cm+^W 

F mn = V«*+( rn-simBlH)* 


/ n £ — "S ~' ' 2??i 

\ B Hi ^ E nn F mn [(m+slB)F mK +(m-slB)E mn \ 

(A 2) 


It is possible to sum the series for o-(0, s/B, BfH) exactly 
only when s/B—^. In this case equations (Al) and (A2) 
yield 


(o4>|)=2Sp ± 

\ — JB./ n~l ft m= — oc 


»*+7 




=2 2 j 

n=I ft ot-+co 


^In*+(m+±y (BIH ) 2 yjn>+(m -±) 2 (B/H )* J 

2 ( m+ l) 

-y/n 2 +(m+^ (B/H) 2 


+2 2 


m 


=*i (m 2 -l/4ty(B/H) 3 


2 “1 
2 




4 ^ 1 4 , 4 11 

hi n 2+ (B/H) 3 ~B/H [_ 6 + (5/S) 2 J 

For other values of s/B it is necessary to sum the series numerically. The series may be rewritten 


(°4!)- 8 si 


m 


2V 

•2S 


m 


E mn F nn [(m+s/B)F mn +(m-s/B)E mn ]^ hi n \l n *+(s/H)^(B/Hf hi [m 2 -(s/Bff 


-R 


.v (0, J- jj) 

(A3) 


where AT is an arbitrary positive integer and l?.v(0, s/B, B/H) 
denotes a remainder term. In evaluating <r(0, s/B, B/H) in 
the present report, N was usually taken equal to 7 and the 
summations indicated in equation (A3) were carried out ex- 


actly; whereas an approximate value was used for the re- 
mainder term. An approximate formula for R.\(0, s/B, B/H) 
is 


R. 


■ w ( 0, i§) 


2^l(B/H) 2 (N+yj+m 


(B/Hf (n+£) (V+f) a |_ 

2^/(B/H) 2 (l/±)+(N+yf 


B/H (B/Hf 


(B/H?(n+ i) 


V(W+i 


(B/Hf(N+±) (B/H) 2 (a 7 +!)J 

It will now be shown how formula (A4) is obtained. It is easily verified that 

L dy f 

J 0 J Q 


dx 


yiB/Hftf+b 2 


and that 


[y 2 +(BIH) 2 xY fl (B/H)b 1 (B/Hf a (B/H)*ab 

dx 1_ 

„ 3?~2a? 


(A4) 


(A5> 

(A6) 
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provided a, 6> 0. Now R N ( 0, slB , B/H) is approximately 
equal to I?, v ( 0, 0 ,B/H) and it is easily found that 


lis(0,0,BIH)=i(± ± + ± S) 

Yn**l m=N+ 1 n=N+l rrt-1.' 


[n i +m t (BlH)T 12 ' 


00 1 9 ” 

2 2 rs+zD/rxis S 

n*AT+l W m=iV-H 


_ 1 _ 

m 3 


These summations may be approximated bj r suitable inte- 
grals so that approximately 


a, (o. g. |)- 4 0," M J/* + J J + , 

2 + 2 ** ' 2 2 

■ L/ v Lj ix ) \ 


n+; “jv+i 


[* 1 (w]J v+ ii* 

■*' 1 • *> 


If the integrals of equation (A7) arc evaluated by means of 
equations (A5) and (AO), then equation (A4) is obtained. 



APPENDIX B 

LIST OF IMPORTANT SYMBOLS 


B tunnel breadth or diameter 

H tunnel height 

C tunnel cross-sectional area 

c chord of airfoil or body of revolution 

t maximum thickness of airfoil or body of revolution 

V volume of wing or body of revolution 

S m maximum cross-sectional area of body of revolu- 

tion 

s half span of wing 

C 0 drag coefficient 

S model area on which drag coefficient is based 

U stream velocity 

M Mach number 

R Reynolds mnnber 

7 ratio of specific heat of gas at constant pressure to 

specific heat at constant volume (c p /c e ) 
p mass density 

q * dynamic pressure 

A, Ki, K% factors depending on shape of airfoil base profile 
(See equations (II), ( 12 ), and (13) and tables 
II and III.) 

X, K 3 , Ki factors depending on shape of body of revolution 
(See equations ( 22 ), (23), and (24) and tables 
IV and V.) 

AO/') linear compressibility correction factor for virtual 
volume (reference 11 ) 

r factor depending on tunnel shape and wing-span- 

to-tunnel-breadth ratio (See equations ( 10 ) and 
( 21 ) and table I.) 

K total blockage correction (See equation (41).) 

K m wing-blockage correction (See equations (42) and 

(43).) 

Ki body-blockage correction (See equations (44) and 

(45).) 

K a k wake-blockage correction (See equation (46).) 
Superscript: 

(') when pertaining to fluid properties, denotes values 

existing in tunnel far upstream from model; 
when pertaining to airfoil characteristics, denotes 
values in tunnel, ^coefficients being referred to 
apparent dynamic pressure q' 


Subscripts: 

(Used only when necessary to avoid ambiguity) 

c denotes values in compressible fluid 

i denotes values in incompressible fluid 

w denotes values for wing 

b denotes values for body of revolution 

wk denotes values for wake. . 

* ■* 

REFERENCES 

1. Lock, C. N. H.: The Interference of a Wind Tunnel on a Sym- 

metrical Body. R. & M. No. 1275, British A. R. C., 1929. 

2. Glauert, H.: Wind Tunnel Interference on Wings, Bodies and Air- 

screws. R. <fc M. No. 1566, British A. R. C., 1933. 

3. Allen, H. Julian, and Vincenti, Walter G.r Wall Interference in a 

Two-Dimensional-Flow Wind Tunnel, With Consideration of 
the Effect of Compressibility. NACA Rep. 781 1944. 

4. Thom, A.: Blockage Corrections and Choking in the R. A. E. 

High Speed Tunnel. R. & XI. No. 2033, British A. R. C., 1943. 

5. Goldstein, S., and Young, A. D.: The Linear Perturbation Theory 

of Compressible Flow, with Applications to Wind-Tunnel Inter- 
ference. R. & XI. No. 1909, British A. R. C., 1943. 

6. Tsien, Hsue-shen. and Lees, Lester: The Glauert-Prandtl Approxi- 

mation for Subsonic Flows of a Compressible Fluid. Jour. 
Aero. Sci. vol. 12, no. 2, April 1945, pp. 173-187 and 202. 

7. von Baranoff, A.: Zur Frage der Kanalkorrektur bei kompressibler 

Unterschallstromung. Deutsche Luftfahrtforschung Forschungs- 
bericht Nr. 1272, 1940. 

8. Gothert, B.: Windkanalkorrekturen bei hohen Unterschalige- 

schwindigkeiten unter besonderer beriieksichtigung des geschlos- 
senen Kreiskaoals. Deutsche Luftfahrtforschung Forsehungs- 
bericht Nr. 1216, 1940. 

9. Gothert, B.: Windkanalkorrekturen bei hohen Unterschallgeseh- 

windigkeiten. Lilienthal-Gesellschaft fur Luftfahrtforschung, 
Bericht 127, Sept. 1940, pp. 114-120. 

10. Wieselsberger, C. : Windkanalkorrekturen bei kompressibler Strom- 

ung. Lilienthal-Gesellschaft fur Luftfahrtforschung, Bericht 
127, Sept. 1940, pp. 3-7.. 

1 1. Ludwieg, H. : Widerstandskorrektur in Hochgeswindigkeitskanalen. 

Deutsche Luftfahrtforschung Forschungsbericht Nr. 1955, 1944. 

12. Gothert, B.: Ebene und raumliche Stromung bei hohen Unter- 

schallgesehwindigkeiten (Erweiterung der Prandtlsehen Regel). 
Lilienthal-Gesellschaft fur Luftfahrtforschung, Bericht 127, 
Sept. 1940, pp. 97-101. 

13. Herriot, John G.: The Linear Perturbation Theory of Axially 

Symmetric Compressible Flow, With Application to the Effect 
of Compressibility on the Pressure Coefficient at the Surface 
of a Body of Revolution. NACA RXI A6H19, 1947. 

14. Young, A. D. and Owen, P. R.: A Simplified Theory for Streamline 

Bodies of Revolution, and Its Application to the Development 
of High Speed Shapes. R. & XL No. 2071, British A. R. C., 1943. 

783 



